Edge states protected by chiral symmetry in disordered photonic graphene 



cn 

o 

<N 

Oh- 

< 

in 



C/5 

o 

• I— I 
-(— > . 

Oh- 

q 
o 

Oh- 



On 
O 

m 



X 



J. M. Zeuner\ M. C. Rechtsman^, S. Nolte\ and A. Szameit^!] 

^Institute of Applied Physics, Friedrich- Schiller- University Jena, Max-Wien-Platz 1, 0774S Jena, Germany and 

^ Solid State Institute and Physics Department, Technion, 32000 Haifa, Israel 

(Dated: April 26, 2013) 

We experimentally investigate the impact of uncorrelated composite and structural disorder in 
photonic graphene. We find that in case of structural disorder not only chiral symmetry, but also 
the vanishing of the density of states at zero energy is preserved. This is in contrast to composite 
disorder, where chiral symmetry as well as the vanishing of the density of states are destroyed. Our 
observations are experimentally proven by exciting edge states at the bearded edge in disordered 
photonic graphene. 



Graphene, a honeycomb patterned monolayer of car- 
bon atoms, has been celebrated as a new wonder ma- 
terial since its first experimental realization in 2004 [1|. 
Especially due to its unique properties such as, e.g., very 
high thermal conductivity, high charge-carrier mobility 
and very high mechanical stability, this material is a 
candidate for many different future electronic and pho- 
tonic devices. Due to its special properties, graphene 
offers a platform for addressing a number of questions of 
a fundamental physical nature. Indeed, novel phenom- 
ena such as the (fractional) quantum Hall effect [2|, |3|, 
strong suppression of quantum interference [^ , and Klein 
tunnelling [5| have been shown to exist in this material. 
The graphene lattice has a two- member basis, composed 
of shifted hexagonal lattices, and can be therefore de- 
scribed in a tight-binding model containing two bands. 
The bands meet at conical intersections (called Dirac 
points) at zero energy, where electrons have a disper- 
sion relation akin to that of the Dirac equation. Impor- 
tantly, the density-of-states (DOS) near those points goes 
to zero. If only nearest-neighbor atoms interact, atoms of 
one sublattice interact merely with the other sublattice 
and not their own. Such systems possess chiral symme- 
try |6l4lO|, from which two important properties follow: 
(1) any state with nonzero energy E must have a chiral 
partner with energy —E and (2) any state with zero en- 
ergy cannot be altered by perturbations respecting the 
chiral symmetry. 

Due to the fundamental difference of Dirac-like disper- 
sion and Schrodinger-like dispersion, a full understanding 
of the impact of disorder on graphene has been recently 
extensively investigated, but a consensus has not been 
fully reached llMl4j . Disorder that is introduced into the 
system may be categorized in two different ways: com- 
posite disorder (also called on-diagonal disorder), which 
corresponds to impurities in the lattice, and structural 
disorder (also called off-diagonal disorder), which cor- 
responds to random displacements of the atoms in the 
lattice. Composite disorder destroys the chirality and 
therefore the Dirac spectrum with the associated vanish- 
ing DOS at zero energy [12|, ll5|, ll6| . Structural disorder 
preserves the chiral symmetry of the honeycomb lattice 
[12|. However, how this affects transport properties and 



the DOS in the Dirac region of the spectrum is still an 
open question. 

In this work, we experimentally probe the impact of 
composite and structural disorder in a graphene-like lat- 
tice, and demonstrate the physical effects for both cases, 
chiral symmetry breaking and preserving. We present 
strong experimental evidence that under structural dis- 
order - which preserves the chiral symmetry - the van- 
ishing DOS at the Dirac point remains. To this end 
we use optical honeycomb waveguide lattices - photonic 
graphene - which constitute an exceptional system for 
probing the unique graphene geometry 171, ll8| . We use 
an optical system since probing the effects of disorder in 
atomic graphene is an essentially impossible task due to 
the difficulty of engineering graphene samples with con- 
trollable amounts of disorder, types thereof, and correla- 
tions therein. In our system, we employ zero-energy edge 
states at the bearded edge of the lattice [19] which turn 
out to be well suited for experimentally analyzing the 
spectrum at the zero-energy region. We find that these 
edge states are destroyed under composite disorder, but 
survive structural disorder. This suggests that the van- 
ishing of the DOS at zero energy is preserved in the lat- 
ter case. Note that such edge states were experimentally 
demonstrated only recently in photonic graphene [20| as 
in electronic graphene an expanded bearded edge was not 
realized so far due to stability problems [21]. A very in- 
teresting feature of such edge states was demonstrated 
lately: In [22, 23] it was shown that edge states expe- 
rience a topological transformation when a linear strain 
is applied to the lattice. In Fig. HJa), a sketch of a 
honeycomb photonic lattice with all basic edge types is 
depicted. A micrograph of the front facet of one fabri- 
cated sample is shown in Fig. [Hb) with the bearded edge 
at the bottom. 

In order to model an optical honeycomb waveguide lat- 
tice, we employ the tight-binding-equations 

<m> 

(1) 

where the summation is taken only over nearest-neighbor 
sites. Here, ifju is the amplitude of the guided mode in the 




FIG. 1. (a) Scheme of the waveguides arranged in a hon- 
eycomb structure, the so-called photonic graphene with the 
labelling of the different edge types, (b) Microscopic image 
of the structure, (c) Dispersion relation of photonic graphene 
which equals that of electronic graphene. 



nth waveguide, z is the propagation direction, Cn^m is the 
coupling constant between the nth and mth waveguide, 
and K, denotes the propagation constant. Evidently, this 
can be written using a tight-binding Hamiltonian matrix 
Htb^ which is identical to that of electronic graphene. 
The only difference of both systems is the evolution coor- 
dinate: whereas in electronic graphene the wave function 
evolves in time, in the optical system the light propaga- 
tion occurs along z. Hence, the temporal dynamics of 
electrons in graphene is emulated by light diffraction in 
the optical waveguide system. In order to achieve time- 
independent equations for the eigenvalues /3, one substi- 
tutes the plane wave solutions tl^{z) = ipe^^^ into Eq. ([1]), 
which results in an eigenvalue equation for the eigen- 
modes of the system: 



■ fill) = HtbiI^ . 



(2) 



For periodic boundary conditions, this equation allows 
the determination of the band structure of a honeycomb 
lattice, where /3 is a function of the transverse Bloch 
wave numbers k = {kx^ky)^ as shown in Fig. [I{c). In- 
deed, there is no bandgap, as both bands intersect at 
six singular points. In the vicinity of these Dirac points, 
the Hamiltonian can be approximated by the relativis- 
tic Dirac Hamiltonian J^d ~ k • <j with a being a vector 
of the Pauli matrices, which is an analogue to the Dirac 
equation for massless relativistic particles [2^. This is 
the central characteristic of the graphene structure and 
one of the main reasons for many exotic properties of 
this material. A further key feature of graphene is the 
observation that at the Dirac points the DOS vanishes. 
In Fig. [2ja) we show all eigenvalues /3 in ascending or- 
der of a system, where periodic boundary conditions are 
applied in all directions. The analytic expression of the 
DOS - which is inversely proportional to the slope of 
the eigenvalue function shown in Fig. [2fa) - is plotted in 
Fig. Efb) and exhibits a clear drop at /3 = 0. The ques- 
tion arises how stable the vanishing of the DOS is under 
the influence of disorder. 



We start our analysis using simulations and exam- 
ine a photonic graphene structure under the influence 
of both different kinds of disorder: composite disorder 
(which manifests in random depth of the lattice sites) 
and structural disorder (which results by a random shift 
of the lattice sites around their average position). Note 
that in our optical system, we can realize both disorders 
independently of each other. Moreover, the structural 
disorder is fully uncorrelated, in contrast to electronic 
graphene, where a random shift of the atoms always in- 
fluences the neighboring sites. We implement the struc- 



tural disorder in Eq. ([I]) by c = cq 



• ^, where cq is 



the mean coupling constant, Ac defines the strength of 
the disorder in units of the coupling constant, and ^ is 
randomly distributed in the sequence [—1,1]. Compos- 
ite disorder is realized by a random distribution of the 
propagation constants: /^ = ^ • ^, where Ak defines the 
strength of the disorder in units of the propagation con- 
stant. Additionally, we assume periodic boundary con- 
ditions in the direction perpendicular to the zigzag and 
bearded edges. In the upper row of Fig. [2] we plot the 
eigenvalues in ascending order for no disorder (Fig. [2fa), 
with a coupling constant cq = 1 cm~^), structural dis- 
order (Fig. Efc), Ac = 8 cm~^^ cq = 5 cm~^ ^ and fixed 
propagation constant kq = 0) and composite disorder 
(Fig. [2fe), An = 8 cm~^ ^ kq = 0, and fixed coupling 
constant cq = 1 cm~^). In the lower row of Fig. [21 we 
plot the according DOS. In order to take into account the 
statistical nature of disorder [25| , in both disorder cases 
we averaged over 5000 realizations, each containing 750 
waveguides. As expected for graphene, the eigenvalue 
spectrum of the completely ordered system shows a DOS 
equal to zero at /3 = (the Dirac points). In case of 
composite disorder the picture completely changes - the 
DOS is not longer zero at /3 = 0. Indeed, with structural 
disorder one still observes a vanishing DOS at /3 = 0. 
Therefore, the chiral symmetry has caused the preserva- 
tion of the Dirac-like DOS even for very strong disorder. 
To further elaborate our results, in Fig. [3] we plot the 
simulated averaged DOS at j3 = for increasing com- 
posite and structural disorder, respectively. The small 
fluctuations are due to the finite number of realizations, 
which were taken for the averaging process. With in- 
creasing composite disorder the DOS at /3 = is increas- 
ing, so that the zero DOS disappears. In contrast, for 
increasing structural disorder the DOS at /3 = stays 
constantly at zero since the chirality is preserved, even 
for very high disorder strengths. We attribute this be- 
havior to the intrinsically different impact of composite 
and structural disorder on the distribution of eigenval- 
ues 26|, l27| . For structural disorder the eigenvalues ap- 



pear in pairs with identical absolute values and opposite 
sign, such that they are symmetrically distributed around 
/3 = 0, which is not the case for composite disorder. This 
is a direct result of chiral symmetry breaking. 

To test the conjectures based on the numerical sim- 
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FIG. 2. Eigenvalues and DOS plotted for three different cases: no disorder, structural and composite disorder. 
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FIG. 3. Dependence of the DOS at /3 = on the strength of 
disorder plotted for structural and composite disorder. 



ulations, we examine the edge states of the honeycomb 
lattice. It is well known that in a finite honeycomb lat- 
tice, edge states exist at the bearded and the zig-zag 
edges (for the nomenclature of the edges see Fig. HJa)). 
Such states exist directly at /3 = 0, and therefore in the 
same region as the Dirac points, where the DOS van- 
ishes. Probing the effects of disorder on the edge states 
of the system represents an ideal diagnostic of the effects 
of chiral symmetry breaking and preservation (e.g., the 
effects on the DOS and the transport), for the following 
reasons. In the case of composite disorder, zero-energy 
edge states may couple to bulk states as a result of the 
disorder, and therefore spread into the bulk. However, 
in the case of structural disorder, the chiral symmetry 
is preserved and the edge states are not perturbed, and 
thus remain closely localized on the edge. Thus, edge 
state properties yield a full description of Dirac dynam- 
ics near /3 = under the influence of disorder. In this 
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FIG. 4. Simulation of the light distribution after z = 7Lc 
when a broad Gaussian beam is launched at the armchair 
edge. The averaging was taken over 200 realizations. The 
edge state in the case of structural disorder propagates in a 
stable fashion, whereas it is destroyed for composite disorder. 



vein, the existence of an edge state can be used as an 
indicator whether the DOS vanishes at /3 = or not. 

The predicted properties of the edge states under the 
influence of disorder are confirmed by tight-binding simu- 
lations of Eq. ([!]), which are shown in Fig. HI The bearded 
edge of the lattice is excited using an elliptical gaussian 



beam. Without disorder the edge state does not spread 
into the lattice. This is because it exists at the same en- 
ergy as the Dirac points where the DOS is zero, so that 
there are no states the edge state could couple to. If 
structural disorder (in the coupling coefficients) is imple- 
mented the edge state is still present and can be observed 
even after a propagation length of z = ILc with the cou- 
pling length Lc = ^ and c as the mean coupling con- 
stant. In contrast, for composite disorder the edge state 
is destroyed and the light spreads into the array as it is 
no longer localized at the edge. Note that in all cases, 
averaging was performed over 200 random realizations. 
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FIG. 5. Experimental light distribution at the end facet of 
the sample, when a broad edge state at the bearded edge is 
excited. The shown images are averaged over 20 realizations. 



In order to experimentally demonstrate the described 
effects, we fabricate a photonic graphene waveguide ar- 
ray using the femtosecond direct- writing technology [28| . 
The structural disorder is realized by randomly shifting 
the positions of the waveguides in the transverse direc- 
tions with respect to their average positions by ±1 jam 
for a disorder of Ac = 0.45 cm~^^ and ±2 jj.m for a disor- 
der of Ac =1.1 c7Ti~^, while the mean nearest-neighbor 
distance was 20 /j.m. The composite disorder was real- 
ized by changing the translation velocity of the writing 
beam during fabrication. The mean velocity for all re- 
alizations was 80 ^^^ with a random shift of ±10 ^^ 
and ±15 ^^ for A/^ = 2.5 cm-^ and /\k = 4.5 cm-\ 

mm ' 

respectively. For the investigation of the light propa- 
gation a broad elliptical beam with its main axis along 
the edge was launched into the bearded edge waveguides 
and the output facet of the sample was imaged onto a 
CCD camera. For each disorder level, we averaged over 



20 reahzations, and the resulting averaged output distri- 
butions are shown in Fig. [5l It is clearly seen that for 
increasing composite disorder the edge state couples to 
the bulk modes and spreads into the lattice. Only for suf- 
ficiently high disorder Anderson localization takes place. 
In contrast, for structural disorder our experiments un- 
equivocally show that the edge state is stable and does 
not couple to the bulk modes. These results represent 
strong experimental evidence that the chiral symmetry 
preserves the vanishing of the DOS at the Dirac point 
{j3 = 0), whereas the breaking of chiral symmetry de- 
stroys it. Note that we checked experimentally that at 
disorder level of Ac = 0.45c?n~^ no Anderson localization 
takes place. 

In conclusion, we theoretically and experimentally an- 
alyzed the impact of composite and structural disorder 
on the DOS of photonic graphene. Our results are gen- 
eral and hold also for electronic graphene, where such 
controlled investigations of disorder are impossible. We 
found that for composite disorder, which breaks the chi- 
ral symmetry of the lattice, the DOS does not vanish at 
/3 = 0. However, for structural disorder we could exper- 
imentally demonstrate that edge states still exist, which 
implies that the DOS vanishes at /3 = even for strong 
disorder, so that chirality prevails. Therefore, all prop- 
erties of the graphene structure that arise from chirality 
and the vanishing DOS at /3 = do not break down for 
structural disorder. Our work may provide further in- 
sight into certain properties of graphene that are strongly 
affected by disorder such as, e.g., the value of the minimal 
conductivity. Furthermore, our photonic system offers a 
unique platform for investigating the influence of prede- 
termined disorder on the electronic transport properties 
of graphene. 
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